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ABSTRACT 

It is widely believed that Type la supernovae (SN la) originate in binary systems where a white dwarf 
accretes material from a companion star until its mass approaches the Chandrasekhar mass and carbon 
is ignited in the white dwarf's core. This scenario predicts that the donor star should survive the 
supernova explosion, providing an opportunity to understand the progenitors of Type la supernovae. In 
this paper we argue t hat rotationis a generic sig nature expected of most non-giant donor stars that 
is easily measurable. ()Ruiz-Lapuente et al.ll2004[) examined stars in the center of the remnant of SN 
1572 (Tycho's SN) and showed evidence that a subgiant star (Star-G by their naming convention) near 
the remnant's centre was the system's donor star. We present high-resolution (R ~ 40000) spectra 
taken with the High Dispersion Spectrograph on Subaru of this candidate donor star and measure 
the star's radial velocity as 79 ± 2 kms~^ with respect to the LSR and put an upper limit on the 
star's rotation of 7.5 kms^^. In addition, by comparing images that were taken in 1970 and 2004, we 
measure the proper motion of Star-G to be = —1.6 ± 2.1masyr~^ and = —2.7 ± 1.6masyr~^. 
We demonstrate that all of the measured properties of Star-G presented in this paper are consistent 
with those of a star in the direction of Tycho's SN that is not associated with the supernova event. 
However, we discuss an unlikely, but still viable scenario for Star-G to be the donor star, and suggest 
further observations that might be able to confirm or refute it. 

Subject headings: astrometry - techniques: spectroscopic - binaries: close - supernovae: general - 
supernova remnants 



Electronic address: I wkerzend® mso . anu.edu . au I 
Electronic address: 'brian@mso.anu.edu.aul 
Electronic address: nomotoQ astron. s . u-tok yo. ac . j p | 
Electronic address: podsiOastro. ox. ac.uk 
Electronic address: anna® astro. as. utexas.edul 
Electronic address: fesenQsnr.dartmouth.edul 
Electronic address: yong@mso.anu.edu.au 

^ Based in part on data collected at Subaru Telescope, which is 
operated by the National Astronomical Observatory of Japan. 

^ Department of Astronomy & Research Center for the Early 
Universe, School of Science, University of Tokyo, Bunkyo-ku, 



1. INTRODUCTION 

Type la supernovae (SNe la) are of broad interest. 
They serve as physically interesting end points of stel- 
lar evolution, are major contributors to galactic chemical 
evolution, and serve as one of astronomy's most powerful 
cosmological tools. 

It is therefore unfortunate that the identity of the pro- 
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genitors of SNe la is still uncertain. For example, without 
knowing the progenitors, the time scales of SNe la en- 
riching the interstellar medium with iron remains highly 
uncertain. But it is the crippling impact on the cos- 
mological application of these objects which is especially 
profound; it is impossible to predict the consequences of 
any cosmological evolution of these objects or even gauge 
the likelihood of such evolution occurring. 

There is broad agreement that the stars which explode 
as SNe la are white dwarfs which have accreted material 
in a binary system until they are near the Chandrasekhar 
mass, then start to ignite carbon explosively, which leads 
to a thermonuclear detonation/deflagration of the star. 
It is the identity of the binary companion that is cur- 
rently completely undetermined. Suggestions fall into 
two general categories (Ibcn 1993): 

• Single degenerate systems in which a white dwarf 
accretes mass from a non-degenerate companion, 
where the companion could be a main-sequence 
star, a subgiant, a red giant, or possibly even a 
subdwarf. 

• Double degenerate systems where two CO white 
dwarfs merge, resulting in a single object with a 
mass above the Chandrasekhar limit. 

The det ection of circum stellar material around 
SN 2006X (jPatat et al.ll2007f) has provided support for 
the single degenerate model in this case, although the 
lack of s ubstantial hydrogen in several other SNe la 
(jLeonardl 12007) poses more of a challenge to this sce- 
nario. 

These models also make different predictions for the 
nature of the system following the explosion. In the dou- 
ble degenerate case, no stellar object remains, but for a 
single white dwarf, the binary companion remains largely 
intact. 

In the single degenerate case, the expected effect of 
the SN on the donor star h as be en investigated by 
[Marietta. Burrows, fc Frvxelll (|2000D . who have calcu- 
lated the impact of a SN la explosion on a variety of 
binary companions. ICanal. Mendez. fc Ruiz-Lapuent j 
()2001f ) have explored many of the observational con- 
sequences of the possible scenarios, and Podsiadlowski 
(2003) has presented models that follow both the pre- 
supernova accretion phase and the post-explosion non- 
equilibrium evolution of the companion star that has 
been strongly perturbed by the impact of the supernova 
shell. To summarize these results, main-sequence and 
subgiant companions lose 10 - 20 % of their envelopes and 
have a resulting space velocity of 180-320 kms~^. Red- 
giant companions lose most of its hydrogen envelope, 
leaving a helium core with a small amount of hydrogen- 
rich envelope material behind, and acquire a space veloc- 
ity of about 10-100 kms-^ iPakmor et all ()2008[ ) have 
used a binary stellar evolution code on a main-sequence 
star and exposed the evolved star to a SN la. Their sim- 
ulations show that even less material is stripped due to 
the compact nature of a star that evolved in a binary. 
We will use their result s wher e applicable. 

iRuiz-Lapuente et all ()2004l . henceforth RP04) have 
identified what might be the donor star to Tycho's SN, a 
SN la which exploded in the Milky Way in 1572. These 
authors presented evidence that this star, Star-G by 



their naming convention, is at a distance consistent with 
the Tycho supernova remnant (henceforth SNR), has a 
significant peculiar radial velocity and proper motion, 
roughly solar abundance, and a surface gravity lower 
than a main-sequence star. However, Star-G is located at 
a significant distance from the inferred center of the rem- 
nant, and any process that has displaced the star must 
preserve the remnant's nearly perfectly circular projected 
shape. During the final stages of refere eing of this paper 
we we re made aware of the article by iHernandcz et all 
(|2009l henceforth GH09), who used Keck HIRES data 
to better constrain Star-G's stellar parameters, and in 
addition, found an enhancement in Nickel abundance, 
r elative to no r mal m etal rich stars. 

Ilhara et al.l (|2007f ) have looked for Fe absorption lines 
from the remnant, using nearby stars as continuum 
sources, with the hope to better constrain the distance of 
these stars to the SNR. With their technique, stars in the 
remnant's center should show strong blue-shifted Fe ab- 
sorption lines, formed by material in the expanding shell 
of Fe-rich material from the SN, moving towards the ob- 
server. Stars in the foreground would show no Fe absorp- 
tion, and background stars both red- and blue-shifted ab- 
sorption. Their study shows that Star-G does not contain 
any significant blue-shifted Fe absorption lines, suggest- 
ing that Star-G is in the remnant's foreground. However, 
these observations and their analysis, while suggestive, 
cannot be considered a conclusive rebuttal of Star-G's 
association with the remnant; this technique requires a 
significant column depth of Fe which is not guaranteed. 
A lack of Fe column depth may be indicated by the fact 
that no stars were found in the vicinity of the remnant 
that showed both blue- and red-shifted absorption lines. 

To further examine the RP04 suggested association of 
Star-G with the SN la progenitor, we have obtained a 
high-resolution spectrum of the star using Subaru and 
its High Dispersion Spectrograph (jNoguchi et al.l ll998). 

We summarize, in section 2, the observational circum- 
stances of the Tycho remnant and any donor star, and 
argue in section 3 that rapid rotation is an important, 
previously unrealised signature in a SN la donor star. In 
section 4 we describe our Subaru observations. Section 
5 covers the analysis of data and the results of this anal- 
ysis. Section 6 compares the relative merit for Star-G 
being the donor star to the Tycho SN or being an un- 
related background star, and in section 7 we summarize 
our findings and motivate future observations. 

2. OBSERVATIONAL CHARACTERISTICS OF THE TYCHO 
REMNANT AND STAR-G 

RP04 have done a thorough job summarizing the rele- 
vant details of the Tycho remnant. The remnant shows 
the characteristics expected of a SN la based on its 
light curve (measured by Tycho Brahe himself), chcm- 
ical abundances, and current X-ray and radio emission 
(jRuiz-Lapuentd l2004[ ) . In figured) we have overlaid ra- 
dio contours'^ on an optical image and have marked the 
position of the stars mentioned in this and RP04's work. 

Although it is not easy to measure the remnant's dis- 
tance precisely, RP04 estimated Tycho's SNR distance to 
be 2.8±0.8kpc, using the ratio of the SN 1006 and Tycho 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Fig. 1. — Radio Contours (VLA Project AM0347) have been 
overlaid (Gooch 1996) on an R-Band Image (NGS-POSS). The 
cutout is an INT image (see text). The stars marked in the figure 
are mentioned in this work and in RP04's work. 

SNR's angular sizes and their re lative ages, and the direct 
distance measure of SN 1006 bv Winkler. Gupta, fc Lo^ 
(2003). iKrause et al.l (|2008.) have recently shown, from 
a spectrum of a light echo associated with the SN1572, 
that this SN was a normal SN la. Using Tycho's ob- 
served light curve, the properties of SN la as standard 
candles, and an extinction value they find a distance to 
the SN of 3.8l]^'5 kpc. Updating their values for the ex- 
tinction values determined in this paper (section 16. ip . 
as well as usi ng an absolute mag nitude for SN la of 
-19.5 ± 0.25 (|Altavilla et all 120041 ) . we find a distance 
of 3.4+};3kpc. In summary, we believe the remnant's 
distance is poorly constrained, but probably between 2 
and 4.5 kpc. RP04 also report the spectroscopic and pho- 
tometric properties for the bright stars near the center of 
the Tycho remnant and find a uniform value of approx- 
imately E{B — V) = 0.6 for stars more distant than 2 
kpc. GH09 have revised the E{B - V) value for Star-G 
to 0.76. 

In addition, for a select list of stars, RP04 provide ra- 
dial velocities and proper motions. For Star-G, RP04 
report a value of Vr = —99 ± 6 kms~^ for the radial ve- 
locity in the Local Standard of Rest (henceforth LSR), 
a proper motion of Hb = —6.1 ± l.Smasyr"^, /i/ = 
-2.6 ± 1.3masyr-\ log.g = 3.5 ± 0.5, and T = 5750K. 
Using HIRES data GH09 have improved the measure- 
ments of Star-G's stellar parameters, finding Vr ~ —80 
kms-i, \ogg = 3.85 ± 0.3, T = 5900 ± lOOK, and 
[Fe/H ] = -0.05 ± 0.09 dex. We note that llhara et all 
(|2007h have classified Star-G as an F8V star (T w 6250 K, 
loggf^ 4.3. IAller et all [l98^, in significant disagreement 
with the RP04 temperature and gravity. We believe the 
GH09 values are based on by far the best data, and for 
the purpose of this paper, we will adopt their values. 

Based on the observations, RP04 asserted that Star- 
G was located at approximately 3 ± 0.5 kpc - consistent 
with the remnant's distance. They note that this star has 
solar metallicity, and therefore its kinematic signature 
was not attributable to being a member of the Galactic 



halo. They further argued that Star-G's radial veloc- 
ity and proper motion were both inconsistent with the 
distance, a simple Galactic rotation model, and the star 
being part of the disk population of the Milky Way. The 
derived physical characteristics of the system were nearly 
identical to what was proposed by Podsiadlowski (2003) 
for a typical SN la donor star emerging from a single 
degenerate system (e.g., U Sco; also see iHachisu et all 
1996; Li & van den Hcuvcl 1997; H achisu et all Il999al : 
iHan fc Pod siadlowskii i2004t iHanll200§) . The revision in 
the stellar parameters by GH09 leads to different dis- 
tance with a larger uncertainty, but by and large, has 
not altered the conclusions above. Taken in total, the 
data provide a rather convincing case for the association 
of Star-G with the Tycho SN. 

3. RAPID ROTATION: A KEY SIGNATURE IN SN lA 
DONOR STARS 

In the single degenerate SN la progenitor channel, mass 
is transferred at a hi gh rate from a secondary sta r onto a 
white dwarf (jNomo to 1982; Nomoto et alll2007t ). These 
high mass-transfer rates require that the secondary star 
overflows its Roche lobe. Due to the strong tidal coupling 
of a Roche-lobe filling donor, the secondary is expected 
to be tidally locked to the orbit (i.e., have the same ro- 
tation period as the orbital period). At the time of the 
SN explosion, the donor star is released from its orbit, 
but will continue with the same space velocity as its for- 
mer orbital velocity and continue to rotate at its tidally 
induced rate. 



i^orb,2= 60 km s"^ 
VaA,2= 90 km 



(Mq) 

Fig. 2. — The expected rotation rate for a donor star as a function 
of its mass at the time of the explosion. The three curves show the 
results for 3 final space velocities of the donor star (similar to those 
suggested by RP04). It is assumed that the white dwarf has a mass 
of 1.4 M0. 



There is a simple relationship between the secondary's 
rotation velocity (worb,2) and its orbital velocity: 



Wrot 



Ml + Ah 
Ml 



■b,2, 



where f{q) is the ratio of the secondary's R oche-lobe 
radiu s to the orbital separation (e.g., given bv lEggletonl 
Il983f ) and q = M1/M2 is the mass ratio of the compo- 
nents at the time of the explosion. Figure [2] shows the 
rotational velocity as a function of secondary mass for 
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several values of Uorb,2 (consistent with RP04s measure- 
ment, and at the low end of values expected for a sub- 
giant star), where we assumed that the exploding white 
dwarf had a mass of 1.4 M©. 

This estimate is strictly speaking an upper limit, as it 
does not take into account the angular-momentum loss 
associated with the stripping of envelope material by the 
supernova and any bloating due to the supernova heat- 
ing. The latter would reduce the rotational velocity to 
first order by a factor equal to the bloating factor (i.e. 
the ratio of the new to the old radius), but the star would 
likely find itself in a state where its radius and tempera- 
ture was atypical of a normal s tar. 

According to the results of M arietta et al.l (|200Cl[ ) , mass 
stripping is not likely to be significant if the compan- 
ion is a main-sequence star or a subgiant. Further- 
more, following bina ry evolution of a main-sequence star, 
iPakmor et al.l ( 20081 ) have shown that even less material 
is stripped. However, if the companion is a giant, it 
would be stripped of most of its envelope. Such a star 
would not show any signs of rapid rotation since the ini- 
tial giant would have been relatively slowly rotating; e.g., 
if one assumes solid-body rotation in the envelope, the 
rotation velocity at ~ 1 Rq will only be ^ 0.5 kms^^ for 
a pre-SN orbital period of 100 d. Moreover, the mate- 
rial at the surface may have expanded from its original 
radius inside the giant, further reducing the rotational 
vel ocity. However, if the stripping is less than estimated 
bv [Marietta et al.l (|2000| ). then it is possible for the sig- 
nature of rotation to persist for a giant, albeit at a much 
lo wer velocity. 

[Marietta et al.l ()2000f ) also showed that due to the in- 
teraction of the SN blast wave with the companion, the 
secondary may receive a moderate kick of up to a few 
lOkms"^, but this kick is generally much lower than 
Worb,2 and therefore does not significantly affect the re- 
sulting space velocity. 

Finally, we note that the observed rotation velocities 
are reduced by a factor sini, where i is the inclination 
angle. However, because the donor star's rotational axis 
can be assumed to be parallel to its orbital axis, a mini- 
mum observed rotation speed can be computed from the 
observed peculiar radial velocity (observed radial veloc- 
ity minus the expected radial velocity of an object at 
that distance and direction) . It is only if the orbital mo- 
tion (and hence final systemic velocity) is solely in the 
plane of the sky, that sini, and therefore, the observed 
rotation, approaches zero. 

4. SUBARU OBSERVATIONS 

To investigate the rotational properties of Star-G, we 
were granted time with the Subaru telescope. Our ob- 
servations of Star-G were taken in service mode on the 
nights of 2005 10 17 and 2005 10 18. 9 spectra were 
taken with th e Hig h Dispersion Spectrograph (HDS, 
iNoguchi et all 11993) with a resolution of R ~ 40000 
(measured using the instrumental broadening of the 
Thorium- Argon arc lines) , an exposure time of 2000 sec- 
onds each (totalling to 5 hours) and a signal to noise 
ratio of about 10 per pixel (measured at 8300 A with 
0.1 Apixer^) . The HDS features two arms, with each 
arm feeding a 2-chip CCD mosaic. The blue arm covers 
6170 A to 7402 A and the red arm 7594 A to 8818 A. An 
OG530 filter was used to block contamination from light 



blueward of our observing window, and data were binned 
by 4 in both the spatial and spectral directions, resulting 
in a pixel size of 0.1 A (at 8000 A) by 0.55". 

Data were pre-processed using tools provided by the 
HDS team and then bias-subtracted. We created a mask 
from bias and flatfielded frames, where we isolated the 
echelle orders and flagged bad pixel regions. The data 
were flatfielded using internal quartz flats, and the 2-D 
images cleaned of cosmic rays (and checked carefully by 
eye to ensure there were no unintended consequences) 
using an algorithm supplied by M. Ashley (private com- 
munication). The spectrum of each echelle order was 
extracted using IRAF'' echelle routines, with wavelength 
calibrations based around low-order fits of a Thorium- 
Argon arc. Wavelength calibration of each extracted 
spectrum was checked against atmospheric O2, and our 
solutions were found to be ac curate in all cases to within 
1 kms~^ (jCaccin et al.lll985l) . Unfortunately, we lacked 
a smooth spectrum standard star for setting the con- 
tinuum, and we resorted to calculating a median of the 
spectra (6 A window) and dividing the spectra through 
this smoothed median. This unusual method was cho- 
sen over the common approach of fitting the spectrum 
with a polynomial, due to the special characteristics of 
this observation (low signal to noise ratio, and a com- 
plex instrumental response). While this does not affect 
the narrow lines our program was targetting, it does af- 
fect broad lines such as the Ha and the Call IR triplet. 
The final step was to combine all spectra and remove any 
remaining cosmic rays (in the ID spectra) by hand. 

5. ANALYSIS AND RESULTS 

5.1. Rotational measurement 

To attain the rotational velocity of the candidate star, 
we measured several unblended an d strong (but n ot sat- 
urated) Fe I lines in the spectrum (Wehrsel ll974l ). Since 
our spectrum only had a combined signal to noise ratio 
of approximately 10, we added the spectra of the lines 
after normalizing them to the same equivalent width. 
As a reference we created three synthetic spectra (one 
broadened only with the instrumental profile, the oth- 
ers with the instrumental profile and Urot sini of 10 and 
15 kms~^ res pectively) with the 2007 version of MOOG 
(|Snedenl[l973l) . using GH09's temperature, gravity and 
metallicity. We use a standard value of /3 = 3/2 for the 
limb darkening although the choice of this value is not 
critical, which we confirmed by checking our results us- 
ing significantly different values of /3. Figure [3] shows 
the comparison between the synthetic spectra of differ- 
ent rotational velocity and the spectrum of Star-G. We 
have scaled the synthetic spectrum using the equivalent 
width. This comparison indicates that the stellar broad- 
ening (rotational, macro turbulence, etc. ) is less than 
broadening due to the instrumental profile of 7.5 kms~^, 
and therefore we adopt 7.5 kms~^ as our upper limit to 
the rotation of the star. If one were to adopt RP04's 
measurements of the peculiar spatial motion, it could be 
concluded that sini is much closer to 1 than (see the 
end of section [3] for further explanation) and thus that 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. — Six observed Fc I line profiles of Star-G are shown 
on the left panel. The right panel shows the combination of these 
line profiles after normalization to the same equivalent width and 
compares them to the spectrum of the Sun, which is convolved with 
3 different values for the rotational broadening kernel. Star-G does 
not show significant rotation, indicating Urot sini < 7.5 kms~^. 



5.2. Radial velocity 

To determine the radial velocity, we used 63 lines to 
measure the shift in wavelength. We find a radial velocity 
in the topocentric (Mauna Kea) frame of reference of 
vtop — —92.7 ± 0.2 kms~^(the error being the standard 
deviation of 63 measurements) . The conversion from the 
topocentric to the Galactic LSR for our observations was 
calculated to be 13.6 kms^^ (IRAF task rvcorrect) using 
the lAU standard of motion. Including the uncertainty 
in the LSR definition, we find a radial velocity in the 
LSR for Star-G of wlsr = -79 ± 2 kms"!. This is in 
significant disagreement with that reported by RP04, but 
agrees with the revised value published by GH09. 

5.3. Astrometry 

RP04 have measured a significant proper motion for 
Star-G of fib = —6.1 ± 1.3 masyr"^, /i; — —2.6 ± 1.3 
masyr"^. Because Star-G is metal rich, and at a dis- 
tance of > 2kpc, this measurement provides one of 
the strongest arguments for Star-G being the donor star 



to Tycho SN. It is almost impossible to account for this 

proper motion, equivalent to a = 58 (^2Tpc) kms^^ or 

3 times the disk's velocity dispersion of = 19 kms^^, 
except through some sort of strong binary star interac- 
tion. 

However, the HST data present an especially difficult 
set of issues in obtaining astrometry free of systematic er- 
rors. For Star-G these issues include the PSF on the first 
epoch WFPC2 image being grossly undersampled, both 
the ACS and WFPC2 focal planes being highly distorted, 
poor and different charge transfer efhciency across the 
two HST images, and that Star-G was, unfortunately, 
located at the edge of one of the WFPC2 chips, making 
it especially difficult to understand the errors associated 
with it. Smaller issues include the small field of overlap 
between the two images, making the measurement sub- 
ject to issues of the correlated motions of stars, especially 
in the /i; direction. 

To cross-check RP04's proper motion of Star-G , we 
have scanned a photographic plate taken in September 
1970 on the Palomar 5 meter, and compared this to an 
Isaac Newton 2.5 m Telescope (INT) CCD archive image 
(INT200408090414934) of the remnant taken in August 
2004. The Palomar plate has an image FWHM of 1.7", 
and the INT image 0.88". While our images have a much 
larger PSF than the HST images, the images have sig- 
nificantly less distortion, are matched over a larger field 
of view with more stars, have fully sampled PSFs, and 
were taken across nearly an 8 times longer time base- 
line. The photographic nature of the first epoch does 
add complications not present in the HST data. The 
non-linear response of photographic plates causes their 
astrometry to have systematic effects as a function of 
brightness (jCannon et al.|[200ll ). especially affecting ob- 
jects near the plate limit, where single grains are largely 
responsible for the detection of an object. 

The position of stars on the INT imag e were matched 
to th e 2MASS point source catalog (jSkrutskie et al.l 
[2006h to get a coordinate transformation (pixel coordi- 
nates to celestial coordinates) using a 3rd-order polyno- 
mial fit with an RMS precision of 40 mas with 180 stars. 
This fit is limited by precision of the 2MASS catalog and 
shows no systematic residuals as a function of magni- 
tude, or position. Using this world coordinate system 
(WCS) transformation, we then derived the positions of 
all stars on the INT image. The coordinates of 60 un- 
crowded stars on the Palomar plate were matched to the 
INT-based catalog, and a 3rd-order polynomial was used 
to transform the Palomar positions to the INT-based po- 
sitions. The fit has an RMS of 65 mas in the direction of 
galactic longitude, and 45mas in the direction of galactic 
latitude. We believe the larger scatter in the direction of 
Galactic longitude is due to the shape of the PSF being 
slightly non-symmetric in the direction of tracking on the 
Palomar plate. This tracking (in RA, which is close to 
the direction of galactic longitude), causes the position 
of stars to depend slightly on their brightness. This ex- 
planation is supported by a small systematic trend in our 
astrometric data in /x;, not seen in /Xf,, as a function of 
rrifi. An alternative explanation is that the trend in /ij 
is caused by the average motion of stars changing due to 
galactic rotation as a function of distance, which is prox- 
ied by m^. We have used the Besangon Galactic model 
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(|Robiii et al.ll2003[ ) to estimate the size of any such ef- 
fect, and find the observed effect is an order of magnitude 
larger than what is expected. The systemic difference be- 
tween assuming either source of the observed effect is less 
than Imasyr"^ in /i;, and has no effect in our /if, mea- 
surement. In our final proper motions, presented in table 
m we remove the systematic trend as a function of ttir 
with a linear function. 

To measure the proper motion of each star, we ex- 
clude each star from the astrometric transformation fit 
so as not to bias its proper motion measurement. Com- 
paring the stellar positions in the 34 year interval we 
find that these 60 stars show an RMS dispersion tr^, = 
2.1masyr~^, a^^ = 1.6masyr~^. For Star-G we mea- 
sure = —1.6 ± 2.1masyr~^, /if, = —2.7 ± 1.6masyr~^; 
this implies that no significant proper motion is detected. 
We do note that this measurement has a similar precision 
to that of RP04, is consistent with no observed motion, 
and is in moderate disagreement with the RP04 measure- 
ment. 
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Fig. 4. — The astrometric motions of 60 stars measured in the 
Tycho SNR center. The measurements have a RMS dispersion 
of 1.6masyr~^. Shown in grey is the proper motion of Star-G 
measured here and by RP04, showing a moderate discrepancy in 
the two measurements. Our measurement is consistent with no 
proper motion. 



In table [T] we present our astrometric measurements of 
all stars listed by RP04 for which we were able to measure 
proper motions. We also give the apparent magnitudes 
in R (partly measured by this work and partly by RP04) 
and the distance from center 9. Due to crowding caused 
by the relatively poor resolution of the first epoch pho- 
tographic plate, several stars are not included that could 
be measured using HST. We include an additional star, 
not cataloged by RP04, which exhibits high proper mo- 
tion. This high proper motion star, which was off the 
WFPC2 images of RP04, we designate HP-l, and has a 
proper motion of /i; ~ 81.3, /if, — —2.7 masyr"^. Due 
to the distance from the remnant's center, (we estimate 



HP-l would have been located 51" from the remnant's 
center in 1572), we doubt this star is connected to the 
Tycho SN, but we include it for the sake of completeness. 

6. DISCUSSION 
6.1. A Background interloper? 

A previously unrecognized property for many progen- 
itor scenarios is the rapid post-explosion rotation of the 
donor (as described in section [3]). The expected rota- 
tion as calculated in Figure [2] is large compared to that 
expected of stars with a spectral type later than F and 
should be easily observable. We have shown Star-G's 
rotation to be less (wiotsini < 7.5 kms~^) than what 
is expected of an associated star if the companion was 
a main-sequence star or subgiant. A red giant scenario 
where the envelope's bloating has significantly decreased 
rotation could be consistent with our observation of Star- 
G, and this will be discussed in section [621 

The primary basis for which RP04 selected Star-G as 
a candidate for the donor star to the Tycho SN was the 
combination of its large peculiar radial velocity and its 
observed proper motio n. In Figure El we use the Be- 
sangon Galactic model (jRobin et al.ll2003l ) to construct 
an expected set of radial velocities for metal-rich stars in 
the direction of SN1572. 
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Fig. 5. — Besangon model for a metal rich {[Fe/H] > -0.2 ) 
Galactic population between and 7 kpc in the direction of Ty- 
cho SNR (1 = 120.1, b=1.4) with a solid angle of 1 square degree. 
The remnant's distance is represented by the black dashed lines (as 
calculated in section [2J. The contours show the radial velocity dis- 
tribution. Our measured radial velocity corrected to LSR and our 
distance are shown, with their respective error ranges, as the black 
rectangle. The distance range calculated by GH09 are indicated 
by the two solid lines. The observed LSR Vr for Star-G is mildly 
unusual for stars at the remnant's distance, and is consistent with 
the bulk of stars behind the remnant. 



Measuring the distance to Star-G is a key discriminant 
in associating the star to the SN explosion. To improve 
the uncertainty of the distance to the star, due both to 
temperature and extinction u ncertainty, we base o ur dis- 
tance on the observed mx (jSkrutskie et al.l [2006) and 
{V - K) color (RP04) We interpolate ATLAS9 m od- 
els without overshoot (jBessell. Castelli. fc Plej|l998l) to 
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TABLE 1 

Proper motions of stars within 45" of the Tycho SNR center. 



a 


5 


Mi 


Mb 
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[hh:mm:ss.ss] 


[dd:mm:ss.ss] 


[mas yr"-"^] 


r 1 1 

[masyr 


[mag] 


[arcsoc] 


Name 


00:25:20.40 


+64:08:12.32 


-0.90 


-0.56 


17.05 


08.9 


c 


00:25:18.29 


+64:08:16.12 


-4.25 


-0.81 


18.80 


10.0 


e 


00:25:17.10 


+64:08:30.99 


-1.82 


1.78 


16.87 


20.3 


f 


00:25:23.58 


+64:08:02.02 


-1.58 


-2.71 


17.83 


31.1 


g 


00:25:15.52 


+64:08:35.44 


1.94 


0.83 


20.28 


31.4 


r 


00:25:15.08 


+64:08:05.95 


-0.67 


1.49 


18.86 


33.3 


j 


00:25:23.89 


+64:08:39.33 


-0.31 


1.08 


19.20 


33.5 


k 


00:25:14.74 


+64:08:28.16 


2.60 


1.46 


17.45 


33.5 


n 


00:25:14.81 


+64:08:34.22 


4.05 


-2.05 


19.35 


35.0 


q 


00:25:13.79 


+64:08:34.50 


2.32 


1.01 


19.90 


41.3 


s 


00:25:14.59 


+64:07:55.10 


-3.94 


2.35 


19.23 


41.7 


t 


00:25:19.25 


+64:07:38.00 


1.75 


-3.43 


16.86 


42.1 


u 


00:25:22.45 


+64:07:32.49 


81.29 


-2.68 


19.81 


48.7 


HP-1 



find a theoretical V — K and absolute magnitude for 
the GH09's values of tem perature and gravity. Using 
a sta ndard extinction law (|Cardelli. Clavton. fc MathisI 
[Mi) {Av = ■i.l2E{B - V) and Ak/Av = 0.109) 
to match the theoretical and observed colors, we find 
Av = 2.58 ± O.OSmag, Ak = 0.28 ± 0.01 mag, and 
e\b -V)^ 0.84 ± 0.05. To better show the uncer- 
tainties, we present our distance moduli scaled to the ob- 
served and derived values of extinction, temperature and 
gravity. The temperature coefficients were determined 
by integrating blackbodies of the appropriate tempera- 
ture with a filter bandpass and fitting a powerlaw to the 
resulting flux. 



{mv - Mv) =12.93 - 3.12(^(5 -V)- 0.84) 



2.5(logg - 3.85) -f 2.5 log 



M 



I Mr. 



4.688 



{jriK - Mk) =12.93 - 0.275{E{B -V)~ 0.84) 



2.5(logg-3.85) + +2.51og 



M 

IMfT 



2.5 log 



V 5900 



1.937 



Assuming a companion mass of 1 Mq we find a (m — 
A/) = 12.93 ± 0.75 mag. This uncertainty is dominated 
by the precision of logg, and equates to a distance of 
D — 3.9 ± 1.6kpc. Star-G, within the errors, is at a 
distance consistent with the remnant. As seen in Figure 
[51 the observed radial velocity of Star-G is consistent 
with a significant fraction of stars in its allowed distance 
range. We also note that if Star-G is indeed associated 
with the SN, that it is likely that Star-G could have a 
mass considerably less than IMq, due to mass transfer 
and subsequent interaction with the SN, although in this 
case, the distance to the star would still be consistent 
w ith SNR dis t ance. 

Ilhara et al.l (120071 ) looked for absorption due to Fe I 
in the remnant's expanding ejecta for 17 stars within 
the Tycho remnant. No such absorption was seen in the 
spectrum of Star-G, potentially placing it in front of the 



remnant. However, the amount of Fe I currently within 
the remnant is uncertain with predicted column densi- 
ties spanni ng several orders of magnitude (0.02 — 8.9 x 
lO^cm-2; iHamilton k Fesen|[l98l lOzaki fc Shigevamal 
[200l . Therefore, we do not believe the lack of significant 
Fe I 3720 absorption in Star-G to be significant. 

In summary, we find that Star-G's radial velocity, dis- 
tance, and stellar parameters are all consistent with an 
unrelated star, but also with it being the donor star. 
There is disagreement in Star-G's measured proper mo- 
tion. The measurements of RP04 are inconsistent with 
normal disk stars at the known distance and strongly 
point to Star-G being associated with the SN, whereas 
the measurements presented here are consistent with a 
normal disk star, unrelated to the SN. In addition, we 
have shown the rotation of Star-G is low (confirmed by 
GH09; Vy-ot < 6.6 kms~^), arguing against association 
with the SN, as does its off center placement in the rem- 
nant. Finally, GH09 have presented evidence that Star- 
G is strongly enhanced in Nickel, an observation that, if 
confirmed, would strongly point to an association of the 
star with the SN. If either the high proper motion, or 
significant Nickel enhancement can be confirmed, then it 
is likely that Star-G is the SN donor star. Otherwise, we 
believe it is much more likely that Star-G is simply an 
interloper. 

6.2. Star-G as the Donor Star to the Tycho SN 

While the case for Star-G's association with the SN 
is not conclusive, it is intriguing, and we believe it is 
worthwhile to look for a consistent solution assuming the 
association is true. While not apriori probable, a self- 
consistent model can be constructed in which Star-G was 
the companion, as we shall discuss now. 

To make such a model work, Star-G has to be a 
stripped giant that presently mimics a G2IV star. At 
the time of the explosion, the star would have been a 
moderately evolved giant (in a binary with an orbital 
period ~ 100 d). The SN eje cta will strip such a giant 
of almost all of its envelope ([Marietta et al.l [2000f ) due 
to its low binding energy; only the most tightly bound 
envelope material outside the core will remain bound. 
Due to the heating by the SN, even this small amount of 
material (perhaps a few xO.Ol Mq) will expand to giant 
dimensions, and the immediate-post-SN companion will 
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have the appearance of a luminous red giant. However, 
because of the low envelope mass, the thermal timescale 
of the envelope is sufficiently short that it can loose most 
of its excess thermal energy in 400 years and now h ave 
the appearance of a G2IV star (jPodsiadlowski 20031). 

A lower mass for Star-G (0.3 — 0.5 Mq) also reduces the 
distance estimate, and makes the observed radial velocity 
more unusual for stars at this distance. The expected 
spatial velocity depends on the pre-SN orbital period and 
should be in th e range of 30—70 km s~^ for a period range 
of 20 - 200 d (jJustham et al.ll2008D . These velocities are 
consistent with the inferred spatial velocity of the object 
relative to the LSR if Star-G is at the distance of the 
remnant, even if no significant proper motion has been 
measured (see Figure [5]) . 

A stripped-giant companion would link the pro- 
genitor to the syrn biotic single-degenerate channel 
(|Hachisu et al.l [T999bf ) for which the symbiotic binaries 
TCrB and RS Oph are w ell studied candidates. In- 
deed, (jJustham et al.l 12008') argued that the ultracool 
low- mass helium white dwarfs (with masses < 0.3 M0) 
that have been identified in recent years are most likely 
the stripped-giant companions that survived SN la explo- 
sions, which could provide some further possible support 
for such a scenario for Star-G. 

If the association is real, Star-G's displacement to the 
SE of the geometric center of the remnant as defined by 
radio and X-ray observations might be interpreted as be- 
ing due to the remnant's interaction with an inhomoge- 
neous ISM. Deep optical images of the remnant do show 
extended diffuse emission along the eastern and north- 
eastern limbs interpret ed as shock precursor emission 
dGh avamian et al.l[2000[ ). This along with an absence of 
detected Balmer-dominated optical emission along the 
whole of the western and southern limbs suggests a den- 
sity gradient of the local interstellar medium with in- 
creasing density towards the NE. An east- west density 
gradient has also been infer red from detailed ra dio ex- 
pansion rate measurements (jRevnoso et al.|[T997l) . Such 
an E-W density gradient could have led to a more rapid 
expansion toward the west giving rise to a small shift in 
the apparent geometric center away from the SE with- 
out creating a highly distorted remnant. However, there 
are problems with this explanation. Deviations from 
spherical symmetry in both radio and X-ray images of 
the remnant are relatively small ( Revnoso et al.l TlQQTt 
ICassam-Chenai' et al.ll2007l ). and the remnant is most ex- 
tended along the eastern and northeastern limbs, just 
where one finds the greatest amount of extended diffuse 
optical emission. Moreover, the remnant's expansion rate 
appears lowest tow ard the northeast (PA = 70 degrees), 
not the southeast (jRevnoso et al.|[l997f) . Although the 
argument that Star-G's SE displacement from the rem- 
nant's current geometric center is a result of an asym- 
metrical expansion is not strong, it remains a possibility. 

The most conclusive way of confirming a stripped-giant 
scenario for Star-G would be an independent, precise 
measurement of the distance to Star-G which in com- 
bination with measurements of the gravity and effective 
temperature would help to constrain Star-G's mass. Un- 
fortunately, such a measurement will most likely have 
to wait for the advent of the GAIA satellite. Alterna- 
tively, one may be able to single out a stripped giant 
from a normal G2IV star through nucleosynthesis signa- 



tures, specifically evidence for CNO-processed material 
(or other nucleosynthetic anomalies). While a normal 
G2IV star is unlikely to show CNO-processed material 
at the surface, a stripped giant is likely to do so. Un- 
fortunately, the data presented here are not of adequate 
quality to explore the detailed properties of Star-G's at- 
mosphere. 

7. OUTLOOK AND FUTURE OBSERVATIONS 

Presently, we believe the evidence for Star-G's associa- 
tion with the Tycho SN is interesting, but not conclusive. 
A possible scenario if Star-G is the donor star, would be 
that of a stripped giant scenario discussed in section 6. 
However, there are still other stars t hat h ave not been 
adequately scrutinized. Ilhara et al.l (j2007l ) have found 
a star (RP04 Star-E) which may contain blueshifted Fe 
I lines, indicating their association with the remnant. 
Unfortunately, the star has nei ther a significant pecu- 
liar radial velocity (jihara eTaLl 12007: RP04) nor a sig- 
nificant peculiar proper motion (RP04 and confirmed by 
our work; see Table [T]). 

High-resolution spectroscopy of each candidate in the 
remnant's center is necessary to precisely determine each 
star's physical parameters. However, the small observed 
velocities of the remaining stars suggest that the donor 
star would have needed to be a giant at the time of ex- 
plosion. Using RP04's observed values, none of the stars 
in the remnant's center appear consistent with what is 
expected of a giant star as the donor star except possi- 
bly for Star-A. We also note that there is an additional 
star present in archived HST images, not cataloged in 
RP04, offset from RP04's star A by 0.5" E and 0.2" N at 
mv = 16.8, {B — V) = 1.0. This star, near the remnant's 
centre, has a color consistent with an F-star (assuming 
that it is behind the bulk of the line of sight reddening) , 
but it will require adaptive optics to obtain its spectrum 
given its proximity to the 13th magnitude Star-A. This 
star could potentially be a non-giant progenitor. 

If future observations are unable to pinpoint a viable 
donor star, other progenitor scenarios will have to be con- 
sidered. These include the double degenerate scenario, 
or a scenario where there is a long time delay between 
the accretion phase of a donor star onto the white dwarf, 
and the ultimate supernova explosion. 
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